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We generalize our analytic approach (Berent M and Parzyriski R 2009 Phys. Rev. A 80 033834) 
to include the effect of strong depopulation of the initial state when considering the multiphoton 
population transfer and scattering of low-frequency light by a 3-level system in the lambda-type 
configuration. We discuss the quality of the approximations made, i.e., the adiabatic and generalized 
rotating-wave approximations. 
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I. INTRODUCTION 

In an analytic approach [1] , we have recently described 
the spectrum of light scattered by a three-level lambda- 
type system under two conditions. One condition was 
that the incident-light frequency was much lower than 
the separation frequencies between the states situated in 
the ends of the arms of lambda configuration (Figure [lj . 
Thus, the transitions in the arms had multiphoton char- 
acter. The other condition was that of weak depopulation 
of the initial state placed at the bottom of one arm (the 
black dot in Figure [l}. As long as the latter condition 
was fulfilled, our analytic approach succeeded in explain- 
ing all details of the spectra found by numerical integra- 
tion of the Schrddinger equation (number of peaks, their 
internal structures, positions, relative heights and sensi- 
tivity to the light-pulse shape and strength) . Taking the 
wavefunction of the system as Ylj=ibj\j), the basis of 
our approach in [T] was a set of two coupled nonlinear 
equations for the ratios r = 62/61 and p = 63/61 of the 
Schrddinger population amplitudes bj. This set, found 



if = {r 2 -l)fi(i) +W2ir - M(t)p, (la) 
ip = (w 3 i + n(t)r)p-M(t)r, (lb) 

with the parameters in it clarified in Figure [T] was an- 
alytically solved for the case of weak depletion of the 
initial state, i.e., under the condition |r|, \p\ <C 1. Under 
this condition, we first dropped the term M{t)p in Equa- 
tion \1'a\ , making this equation a quadratically nonlin- 
ear Riccati-type equation known from the two-level prob- 
lem [3]. Next, the Riccati equation for r was solved 
in the zero-order approximation by neglecting the term 
quadratic in r. A better solution to the Riccati equation 
for r, describing approximately the effect of nonlinear- 
ity in r, was then found in a perturbative way and this 
solution was substituted to Equation ( lb ) . At a given 
r, Equation ( lb ) is a nonhomogeneous linear first-order 



differential equation with known formal solution [T]. 
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Now we extend our approach [T] to cover the case of 
significant depopulation of the initial state, particularly 
the spectacular case when the population is periodically 
exchanged between the two lower states lying in the bot- 
toms of lambda, leaving the highest state practically un- 
populated. Such a population exchange between states 1 
and 3 was, in fact, obtained in the numerical calculations 
reported by both us (Figure 2(c) in pQ and the corre- 
sponding text) and Caldara and Fiordilino (Figure 2(b) 
in [1]). Specifically, an effective exchange between the 
mentioned states occurred when these states were sepa- 
rated by nearly the energy of two photons. This case of 
the two-photon separation of the lower states is the very 
one we focus our attention on throughout this paper. In 
terms of r and p, the case to be considered means that r 
is still small (\r\ <C 1) but p can be large. 




Figure 1: The system under study with the separation fre- 
quencies UJ21 and 0J23 much higher than the laser frequency 
ujq. Q(t) and M(t) are the instantaneous Rabi frequencies for 
the electric-dipole 1-2 and 2-3 couplings, respectively. 



II. THE PROBLEM IN TERMS OF THE 
RICCATI EQUATION 

In Figure [T] we sketch the system under study. The 
initial state of the system has its cigenfrequency lj\ = 0. 
The parity +(— ) of the 1st state differs from that of the 
2nd state and is the same as that of the 3rd state so the 
only non-zero dipole matrix elements in the system are 



2 



A 1 12 = (l|ez|2) and ^23 = (2|ez|3). Instead of applying 
the variables r and p as in pQ, we now prefer to work 
with the equations for the population amplitudes bj Q] : 



ib x = -fi(t)&2, 

ib 2 = u 2 \b 2 - n(t)&i - 

i& 3 = o; 31 & 3 - M(t)b 2 ■ 



M(t)b 3 



(2a) 
(2b) 
(2c) 



where and M(t) are the instantaneous Rabi frequen- 
cies for the electric-dipole 1 — > 2 and 2—^3 couplings, re- 
spectively, and CJ21 = ^2, ^31 = W3 are the eigenfrequen- 
cies of the other states. Precisely, f)(f) = fi^/(i) cos(w t) 
and M(t) = Mxf(t)cos(u)Qt), where /(i) is a slowly 
varying pulse shape function and cos(woi) refers to the 
laser-field oscillations as fast as the optical frequency 
cj < w 2 i, w 23 . Moreover, Q R = an d M fi = 

are the familiar time-independent one-photon Rabi fre- 
quencies with Eq being the peak electric-field strength. 

Under the conditions W21 uo and \b 2 \ 2 <C 1, we 
shall first reduce Equations ([2]) to an effective two-state 
set. Since the numerical calculations [TJ 0] pointed to 
negligible population of level 2 when one meets nearly 
two-photon resonance between states 1 and 3, we approx- 
imately eliminate b 2 from the set (Equations ([2])). First, 
we formally integrate Equation ( |2b[ ) in the range (0,£). 
Since this equation alone is a nonhomogeneous linear in 
62 first-order differential equation, its formal solution is 

PQ 



iU2lt / (fi(i>i(i') + M(i')6 3 (i'))e^ 2lt 'di'.(3) 



We know that O(0) = M(0) = and that (since w 2 i > 
wo) these instantaneous Rabi frequencies oscillate slowly 
in time with respect to the function e 1 " 21 '. However, we 
a priori assume that the amplitudes b\ and 63 also change 
slowly in time as compared to e lU2lt . Thus, the solutions 
for b\ and 63, if found under this adiabatic assumption, 
will have a limited range of applicability. This range will 
be defined later on. With this adiabatic assumption, the 
leading term in Equation ^ is 



W21 



(n(*)6i(t)+M(t)6 3 (t)) 



(4) 



This 62 (t), when substituted to Equations (2a) and (2c I, 
gives a pair of equations for 61 and 63: 



ib\ 



U21 



n(t)M(t) 



^21 



(5a) 
(5b) 



where 0(t) = f*(n 2 (t)/io 21 )dt and tf>(t) = J *(w 31 - 
(M 2 (t )/uj 2 i))dt . Above, 6 and 4> — W31 represent the in- 
stantaneous Stark shifts of states 1 and 3 caused by direct 
coupling of these states to state 2, while Q.(t)M{t)/oj 2 i 
looks like an effective time-dependent coupling between 



states 1 and 3 through the intermediate state 2. Sub- 
stituting bi — u exp(i9(t)) and 63 = v exp(— i<fi(t)), we 
transform set ([5| into 



in 
1 v 



S(t)v, 
-S*(t)u, 



where 



with the phase 



S(t) 



n(t)M(t) 



-ia(t) 



W21 



(6a) 
(6b) 



(7) 



a(t) = 6(t)+cj>{t) 

= w 3 it + A 4 + A 



f 2 (t )cos(2w i )dt 



w 81 * + A t + —f{t) sin(2w t) 



(8) 



depending on the optical field strength through the com- 
bined static Stark shift A i/r 
rt 



Here, At = 



2lo 21 ■ ~' 

A J Q pit )dt . When writing the final (approximate) 
form for the phase a(t) we have applied the condition 
that the pulse shape function f(t) changes slowly with 
respect to cos(2a;oi). Introducing a new variable x = ^, 
set Q is then reduced to a single equation for this vari- 
able, namely the Riccati-type equation 



a = s(t) x 2 - s*(t), 

= {x 2 -l)S(t)+2ilm(S(t)). 



(9) 



This equation has structurally the same form (but dif- 
ferent S(t)) as the one met in the problem of a 2-level 
system [51 [3] . The corresponding two-level system equa- 
tion reads iR = S(t)R 2 - S*(t), where R = (& 2 /&i)e iW2lt 
and S(t) = Vl(t)e~ lU2lt . In terms of the present x and 
S{t) one obtains from Equations (|6| the required ampli- 
tudes: 



&i (t) = exp 
bs(t) = x exp 



S(t )x{t )dt + 6{t) 



1 I / S(t )x{t )dt - 4>(t) 



, (10a) 
(10b) 



Finally, using the Fourier-Bessel expansion to the term 
e ia ^\ i.e., exp(izsin^) = Jni z ) exp(m^) and 

also the relation J„_i(z) + J n +i(z) — ^J n {z) satisfied 
by the Bessel functions J n (z) we write the present S(t) 
defined by Equation Q as 



S(t) = 



QrMr 
2u 21 



fit) £ (-i)' 



(ii) 



00 

= ^^ e -i[(w3i-2nwo)t+A t ] 
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This S(t) was derived taking fast oscillations of the elec- 
tric field in the optical pulse in the form of cos(cl>o£). For 
the sin(w t) oscillations one should replace (—1)™ by 1 
in Equation (111. Equations §f§ and (11) for S(t) show 
a great role played by the term f 2 (t)sm(2(jJot) in the 
phase a(t). First of all, it is a term allowing higher-order 
even-photon 1 — > 3 resonances in the system, i.e., higher 
than two-photon resonance. Moreover, this term intro- 
duces a complex field dependence to a fixed resonance 
because, due to this term, the field appears in the argu- 
ment of the Bessel functions. In Section III, more details 
will be given on the role of this term in the phase a(t). 

The key Riccati Equation ^ is seen to play a cru- 
cial role in the above presentation. Obviously, it also 
results from Equations ([I]) for the variables r and p. To 
show this, we first neglect the term quadratic in r in 
Equation (la I when considering the case of weak pop- 
ulation of state 2. Then, the simplified Equation ( la 



is solved in the adiabatic approximation with the result 
r ~ (f2(£) + M(t)p)/oj 2 i analogous to Equation Q. Sub- 



stituting this r to Equation ( lb) and taking into account 
that p = le - ™' 1 ' fives the Riccati Emia.tion (191) 



gives the Riccati Equation Q9p . 



where r = u^t. Substituting the above x(t) to Equa- 
tions ( 10 ) results in 



(15a) 



/-,(/) = <■<>* ( I S(t')dtj e l9(t) (3(t). 
/-:!(/) = <J / sinf I S(t)dt 

-I r, cos ( I S{t)dt ) \ < )'(/). (l.-)l>) 



In the following, 



where f3(t) = exp i f Q S(t )x\(t )dt 

we shall consider a special form of these amplitudes. 

Let us assume that the pulse envelope function f(t) 
rises from to 1 in a short few-cycle time t r and then 
keeps constant value f(t) = 1 for t > t r . Thus, for 
the rising time much shorter than the pulse length, we 
roughly approximate A t ~ At. As a consequence, when 
the dynamic even-photon 1 — > 3 resonance takes place 
(w 3 i + A = Pujq, where P = 2n = 0, 2, 4, . . . if w 31 > 0; 
P is negative if W31 < 0), S(t) of Equation (11) splits 
into its slowly and rapidly varying in time parts. The 
slow resonance part is that with n = P/2 only, 



III. APPROXIMATE SOLUTION OF THE 
RICCATI EQUATION 

In [TH3], an analytic solution to the Riccati equation 
like Equation ([9| was given but for the case of weak de- 
population of the initial state 1. That solution is not 
valid for the present case because now the initial-state 
population can completely be moved to the third state 
under the condition of two-photon 1 — > 3 resonance. The 
strong pumping of state 3 means intermediate and large 
values of x 2 and, thus, the dominant role of the first term 
on the right-hand side of Equation On the other 

hand, the second term on the right-hand side of Equa- 
tion (|9| is expected to be non-negligible in the opposite 
limit of weak pumping of state 3, i.e., when x 2 is small. 
Mathematically, if 



x (t) = i tan / S(t )dt 



(12) 



is the solution of the reduced Equation (JoJ) , ix Q = (x 2 
l)S(t), and 



xi{t) = 2 / Im(5(i ))dt 



(13) 



solves the other reduced Equation ([9]), ix\ = 2ilm(S(t)), 
then 



x{t) — x (t) + Xi(t) 



(14) 



is an approximate solution of the complete Riccati Equa- 
tion |9]) under the condition \(x 2 +2xoXi)(S(t)/u>o)\ <§C 1, 



Sr(t) = (-1) 



^ R M Rfm 



2w 2 i 



P 2w 
2A/1 



A/ 2 W 
2w 



(16) 



while the rest in S(t) is the other part. This slow part 
will give dominant contribution to the integral f Q S(t )dt 
rooted as the argument in the trigonometric functions in 
Equations (15). Thus, in analogy to the standard one- 



photon rotating- wave approximation (RWA) [5], we now 
make generalized rotating- wave approximation (GRWA 
[BrE]) in the case of a given even-photon 1 — > 3 reso- 
nance, i.e, neglect the rapidly oscillating terms in S(t). In 
GRWA, S(t) = fl^(t) and is thus real entailing x x = 0, 

= 1 and x = x = itan ( f* fi( p )(i')dt'^ . The last x 

is seen to be the exact solution of the Riccati Equation ^ 
in its GRWA limit. The GRWA population amplitudes, 
obtained from Equations (15), are thus 



b^t) = cob (J n (p \t)dt^j e l9{t) , (17a) 

bs{t) = * sin (^ V {P) (t)dA e _< *W. (17b) 

After the rising time t r has passed, the GRWA popula- 
tions of levels 1 and 3 (|&j| 2 ) are seen to oscillate at a 
constant P-photon Rabi frequency resulting from Equa- 
tion (16) by putting /(<) = 1. 

The applicability range of the above special solution 
(Equations ( |17[ )) results from the requirement that 61 and 
&3 change slowly in time with respect to exp(iw2i£), where 
u 2 i ujq. The amplitudes will change slowly if four 
inequalities are fulfilled: |w 3 i/cj 2 i| -C 1, (£Ir/uj2i) 2 -C 1, 



4 



(M R /uj2i) 2 < 1 and \il {p ^ /u 21 \ < 1. The first inequality 
is realized in the lambda-configuration of the levels only 
(Figure [I]), while all other inequalities impose the upper 
limit on the applied electric field. Since W21 3> wo, the 
strongest field is such that both Q r /ujq and Mr/ojq are 
much smaller than ^21/^0 allowing Qr/ljq and M r /loq to 
be as high as of the order of 1. For this strongest field, one 
finds |A/2a>o| <C 1 justifying the approximation for the 

(A/ 3 (t)/^ ) P /a 



Bessel function Jp/ 2 ( ^ 



(P/2) 



With 



this approximation, Equation (16) for fl^ p \t) reduces to: 



photon order P. For P = 0, i.e., the case of degen- 
eracy (w3i = Sp — > 0) there is no intermediate state 
between the initial and final states characterized by ceo 
and /3q, respectively. There is also no intermediate state 
for P = 2 but now the initial and final states are 
characterized by ctQ and /3i, respectively. However, for 
P = 2m > 4 there are P — 2 intermediate states char- 
acterized by ax, Pi, a 2 , Pi,..., a(p/2)-i, P(p/2)-i- By 
adiabatic elimination of all intermediate states one ob- 
tains a pair of equations for the amplitudes of the initial 
and final Floquet states coupled by P-photon resonance: 



n(0) = ^ /2(t); (18a) 
2o; 2 i 

( w y 

n (p>o )(t) = ( _ 1) , n w (t) i_^L i (18b) 



where q = ^ — 1. 

Nearly the same solution as our Equations ( fl7| ), but 
with rt p \t) given by Equations (I18J) rather than gen- 
eral Equation ( [16] ), can be obtained from the Floquet-like 
analysis. However, this analysis seems to be less comfort- 
able and consumes much more time than our derivation. 
The basis for this analysis is the series expansion for the 
atomic amplitudes bx (t) and 63 (t): 

00 

bx(t) = a m {t)e- am ^\ (19a) 

m. — — 00 
00 

b 3 (t) = J2 Pm(t)e-* 2muJot . (19b) 

m— — 00 

When substituted to Equations d5]), this expansion re- 
sults in the evolutions for the Floquet-state amplitudes: 



!«0 



2uj 2 



P/2 



= \5p- M «« 



2w 2 i 
rt p \t)a , 



fit) Pp/ 



(21a) 



(21b) 



where ^ p \t) is given by Equations ([18]). The last set 
looks structurally as Equations |5| and, under the con- 
ditions that led us to Equations (17), we have 



bx(t) 



hit) 



cos 

exp 

/3p/2e 



exp 



3l 

2^21 Jo 



fit )dt 



(22a) 



n (p \t)dt 



i[At - Puj t 



.(22b) 



n 2 

-(2muj + —^f(t))a m 

2UJ 2 1 

tt 2 

T^f(t)ia m -i +a m +i) 

^M^f(t)(2(3 m + p m _x + p m+1 ) , 
4a; 2 i 



(20a) 



(oj 3 x - 2mw 



M 



R f 2 



2u- 



21 



fitWn 



M 2 

Y^fWrn-l+Pm+l) 

4w 2 i 



(20b) 



QrMr 2 



4w 2 



fmc 



a m -x + ot m+ x) ■ 



Under the condition of P = 2m-photon resonance (W31 — 
Pljq = 5p — > 0), Equations ( [20] ) can be solved approx- 
imately in a truncated basis of states. The states of 
this basis are those forming the shortest coupling chain 
between the initial state of amplitude ckq and the final 
state of amplitude P m —p/ 2 . The number of intermedi- 
ate states in this shortest chain depends on the multi- 



The advantage of our Equations (17) is that they in 



elude the Rabi frequencies in their general form given by 
Equation (16), while the solution from the Floquet-like 



analysis (Equations (22)) includes the approximate Rabi 



frequencies given by Equations (18). 



IV. RESULTS 

In Figure [2] we show the quality of the analytic GRWA 
population amplitudes (17). We compare the results 
obtained from Equations (17) (solid lines) with those of 
the numerical integration (dotted lines) of the starting 
Equations @. To ensure roughly the applicability 
of the analytic solution and the 1 — > 3 even-photon 
resonance we have taken the same parameters as in 
DP, i-e., w 2 i/a;o = 13, W31M) = 1.99445, fl R /u) = 0.5, 
Mr/wq = 0.3 and, additionally, 4-cycle rise time of the 



pulse (/ (t) = sin 



g- ) for < t < t r = 4T = ^ and 

f(t) = 1 for t > t r = 4T). The assumed u 2 x, Q R , and 
M R give A/wo = 0.00615 so one meets the dynamic 
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Figure 2: (Color online) (a) Evolution of the real and imag- 
inary parts of the initial-state population amplitude bi (we 
do not show the evolution of 63 due to fast oscillations with 
frequency 2tJo under slowly changing envelope), (b) Two- 
photon Rabi oscillations of the populations in states 1 and 
3. (c) Coherent spectrum of scattered light. The conditions: 
= 13, wsi/wo = 1.99445, n R /uj = 0.5, Mr/uo = 0.3, 
the pulse shape function f(t) increases to its maximum value 
of 1 by four optical cycles and then keeps this value. Solid 
lines - the analytic results from Equations (171, dotted lines 



the results of numerical integration of Equations |2| 



2-photon 1 — > 3 resonance (W31 + A = 2.0006u;o). 
For the above frequency-strength parameters the 
analytic amplitudes ( |17| are reduced to the forms 

h ~ cos(fi( 2 ) f*f 2 (t')dt) ■ exp(^ f*f(t')dt 
and 63 - ismfnW J*f(t')dA 



Figure 3: (Color online) (a) Evolution of the populations in 
states 1 and 3. (b) Coherent spectrum of scattered light. The 
conditions as in [4]: W21/W0 = 19, 0)23/1^0 = 17, lo-ji/ujo = 2, 
Q r /lj = 0.8, Mr/ujo = 0.7, f(t) = ^-f exp [l - (^) 2 ] with 

t > and t = 100 T. (c) As (a) but now for W21/W0 = 
50, u)2z/ojo = 48 and the same all other parameters. Solid 
lines - the results after the adiabatic and generalized rotating- 
wave approximations, dotted lines - the result of numerical 
integration of Equations |2| . 



exp 



where ^ = 



results from the approximate Equation |l8b| for Q( 2 \t) 
We have used these forms of b\ and 63 when drawing 
all analytic curves (solid lines) in Figure [2] Figure [2^ a) 
shows the evolutions of the real and imaginary parts 
of b\, Figure f2Vb) the 2-photon Rabi oscillations of 
populations (|&i] 2 ,|&3| 2 ) between levels 1 and 3 (the 
population of level 2 does not exceed 0.0015 for the 



G 



parameters chosen), and Figure [2^c) the coherent 
part of the spectrum of light scattered by the sys- 
tem (S c (uj) means squared modulus of the finite-time 
Fourier transform of the induced dipole in the system: 
dinduced = 2/ii2-Re(6*6 2 ) + 2^23^(6263), where 6 2 is to 
be calculated from Equation Q by applying the above 61 
and 63). As seen, the analytic solution explains well the 
trends observed in the results of numerical integration 
of Equations ph. Moreover, much better coincidence is 
now found between the analytic and numerical spectra 
in the low-frequency regime than in our previous paper 
(Figure 2(c) in pQ). This low-frequency part of spectrum 
includes three components. The strongest component 
simply corresponds to the elastic Rayleigh scattering. 
The two remaining components, 3wo and 5wo, are the 
result of nonlinearity in the interaction between the 
system and light. As the field intensity increases (large 
r, Qr and A/r), the 5wo component becomes more 
pronounced, as seen in Figure |3jb). One would expect 
the appearance of higher-order odd harmonics of ujo in 
the spectrum as the intensity rises. However, we would 
like to mention that, in the high-frequency regime of 
the numerical spectrum in Figure [2](c) (a; ~ W21), one 
observes an additional low structure. The height of this 
structure is about five orders of magnitude lower than 
that of the leading peak in Figure [2]Jc). Our analytic 
solution was not able to reproduce this additional low 
structure observed in the numerical spectrum. This 
is, probably, the cost we pay for the approximate 
procedures, particularly the procedure of adiabatic 
elimination of 62 from Equations ([2]). The adiabatic 
elimination has turned out to be the main reason of 
quantitative differences between the solid and dotted 
curves in Figure [2j It was found that this procedure 
works better for higher ratios u>2i/^o- By increasing this 
ratio well above the assumed value 13, the differences in 
the long-time population evolutions gradually decrease. 
Obviously, longer pulses are needed to observe the 
2-photon Rabi oscillations when increasing the w 2 i/wo 
ratio. 

In Figure |3j we present a different comparison for the 
conditions assumed by Caldara and Fiordilino |3] in their 
numerical integration of the Schrodinger equation. These 
conditions were: W21/W0 = 19, ujsi/luo = 2, VLr/ujq ~ 0.8, 
Mfi/u>o = 0.7 and the smooth pulse shape function 

/(t) = (A ) 2 exp[l-(^ 



where t > and r = 100 T. 



ically the GRWA version of Equations ^ for u and v 
(S(t) — > QS p \t)) and then obtained the solid-line curves 
shown in Figure [3] These curves, resulting from both 
the adiabatic and generalized rotating-wave approxima- 
tions, are compared with the dot-line curves coming from 
the numerical integration of Equations ([2|. The com- 
parison shows that, under the assumed conditions, the 
two approximations lead to results reflecting the main as- 
pects of the exact (numerical) solution of Equations ([2| . 
Again, by increasing the ratio W21/W0 (from 19 to 50) 
much better agreement was found between the analytic 
and numerical curves (see Figure [3]jc)). In the same time 
scale as in Figure |3j a), almost all the population is seen 
in Figure [3]jc) to move from the initial state 1 to state 
3 but without the return 2-photon transition from 3 to 
1. The residual differences, still seen between the solid 
and dotted curves in Figure |3jc) , may suggest that the 
approximate GRWA multiphoton Rabi frequencies given 
by Equations ( 18 ) are not sufficiently accurate. rlow- 



For these conditions, Equations ( 16 1 and ( 18 ) for (t) 
need to be multiplied by the factor e~' lAt . After including 
this factor we were not able to solve the Riccati Equa- 
tion ([9]) in an exact analytic way. So, we solved numer- 



ever, taking il^(t) in the form of Equation (16) has 
not removed the residual differences. As to the photon- 
emission spectrum shown in Figure [3][b), we do not ob- 
serve now any low structure in the high-frequency regime 
(uj ~ W21) of the numerical spectrum (dotted line). This 
makes a difference when compared to the numerical spec- 
trum in Figure ^c). 



V. SUMMARY 

With this paper we have extended the scope of valid- 
ity of our recent analytic approach [T] to the problems of 
multiphoton population transfer in a three-level lambda- 
type system and scattering of low-frequency light by this 
system. The previous approach covered the case of weak 
depopulation of the initial state only. Now, the present 
version of the approach allows strong depopulation of the 
initial state as well. The amended approach was shown to 
explain the main numerical results, particularly the two- 
photon Rabi oscillations of population between the states 
forming the bottoms of A and the dominant (if not all) 
peaks in the spectrum of scattered light. Though we con- 
sidered the 2-photon resonance between the lower states 
in A, the approach holds for any even-photon resonance 
between the mentioned states as well as for the case of 
degeneration of these states. 
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